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EXECUTIVE  SUMMARY 


Space-based,  synthetic  aperture  radars  (SAR)  are  being  considered 
for  a  variety  of  strategic  and  tactical  military  missions.  SARs  in  space 
can  achieve  resolutions  comparable  to  practical  opticaL  instruments,  and 
they  are  unaffected  by  clouds,  weather,  and  darkness.  The  major  obstacle 
to  their  use  is  the  extensive  data  reduction  required,  although  this 
problem  is  common  to  any  advanced  surveillance  system  and  rapid  progress 
in  overcoming  it  is  being  made. 

The  azimuthal  resolution  of  a  SAR  with  a  fixed  antenna  is  limited 
by  the  beamwidth.  Practical  space-based  systems  can  synthesize  apertures 
approaching  %  20  km.  At  L-band ,  the  corresponding  angular  resolution 
is  a  few  hundreths  of  mi 1 l irad ians ;  however,  naturally  occurring  iono¬ 
spheric  disturbances  can  cause  comparable  amounts  of  angular  jitter. 

Viewed  another  way,  the  spatial  coherence  scale  of  an  L-band  wave  passing 
through  a  disturbed  ionosphere  can  be  reduced  below  20  km. 

To  investigate  the  potential  impact  of  the  disturbances  on  SAR 
imagery,  we  have  applied  a  general  space-based  radar  model  developed  in 
the  first  phase  of  this  contract.  The  resolution  characteristics  of  a 
SAR  are  determined  by  the  range-azimuth  ambiguity  function.  By  simulating 
ionospheric  disturbances,  we  showed  that  moderate,  naturally  occurring 
disturbances  degrade  the  azimuthal  resolution  primarily  by  raising  the 
s ide lobe  le  ve  1  . 

To  look  for  evidence  of  the  potential  effects  of  such  disturbances, 
we  reviewed  SAR  data  from  the  SEASAT-A  satellite.  One  pass,  which  was 
coincident  with  incoherent-scatter  radar  operations  that  measured  a 
highly  structured  F-region,  showed  an  image  region  with  a  conspicuous 
contrast  reduction.  Other  evidence  of  such  effects  during  disturbed 
conditions  also  have  been  reported. 
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We  believe  these  results  to  be  the  first  demonstrated  evidence 
of  propagation-induced  SAR  image  degradation.  The  results  are  consistent 


with  predictions  based  on  the  SATCOM  channel  model  that  has  been  used 
extensively  for  analyzing  communication  and  navagation  systems.  The 
inputs,  however,  draw  heavily  on  our  structure  phenomenology  program 
for  the  natural  ionosphere. 

More  detailed  simulations  are  planned  to  identify  the  characteristics 
of  severe  disturbances.  This  final  report  completes  the  development  of 
the  theoretical  model. 
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INTRODUCTION 
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Space-based  radar  (SBR)  systems  have  been  proposed  for  future  CONUS 
defense,  surveillance,  and  tactical  battle-support  functions.  These  systems, 
however,  may  have  to  operate  in  naturally  or  nuc 1  ear -d  i  s tur bed  propagation 
environments.  This  final  report  describes  Lhe  effects  of  propagation 
disturbances  on  space-based  synthetic  aperture  radars  ( SARs ) ,  which  are 
becoming  increasingly  more  important  for  strategic  and  tactical  support 
m i ss ions  . 

In  Topical  Report  L,^  we  developed  a  general  mathematical  formalism 
for  using  the  satellite  communications  model  to  simulate  SBR  propagation 
effects  and/or  to  calculate  measures  of  average  performance  degradation. 

The  functional  elements  of  a  monostatic  radar  system  are  diagramed  in 
Figure  l,  which  is  taken  from  Topical  Report  1. 


TARGET 


FIGURE  1  FUNCTIONAL  DIAGRAM  OF  A  MONOSTATIC  RADAR  SYSTEM 


References  are  listed  at  the  end  ot 


the  report. 
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The  SATCOM  model  characterizes  the  propagation  effects  on  a  single 
path  in  terms  of  a  time-varying  impulse  response  or  the  equivalent  trans¬ 
fer  function.  To  accommodate  a  two-way  path  fully,  the  SATCOM  channel 
model  effectively  must  be  evaluated  twice,  although  the  differences  between 
the  paths  can  be  accommodated  by  scaling  the  independent  variable  appro¬ 
priately. 

The  computation  of  average  performance  measures  for  a  two-way  path 
can  be  very  complicated.  However,  in  Topical  Report  1  we  showed  that, 
under  the  severely  disturbed  conditions  of  primary  concern,  the  forward 
and  return  paths  are  uncorrelated.  Thus,  average  performance  measures 
generally  can  be  constructed  from  products  of  average  measures  on  the 
individual  paths. 

In  addition  to  the  two-way  path  problem,  propagation  effects  in 
SBR  systems  are  more  critically  dependent  on  angular  deviation  than  are 
propagation  effects  for  SATCOM,  where  the  source  (target)  location  generally 
is  known.  This  is  particularly  true  for  SARs. 

The  model  developed  in  Topical  Report  1  uses  a  discrete  array  of 
noninteracting  point  sources  to  simulate  an  arbitrary  antenna  configuration. 
For  a  dense  distribution  of  radiators  (or  receivers),  the  more  familiar 
aperture  distribution  functions  can  be  recovered  by  replacing  summations 
with  integrals. 

The  actual  propagation  effects  are  accommodated  by  applying  the 
appropriate  channel  transfer  function  to  the  radar  signal.  For  a  narrow- 
band  signal,  the  propagation  effects  are  multiplicative.  At  any  instant 
in  time,  the  effect  is  to  convolve  the  gain  pattern  of  the  antenna  with  a 
random  function.  The  average  distortion  of  the  gain  pattern  is  obtained 
by  computing  the  Fourier  transform  of  the  aperture  distribution  multiplied 
by  the  mutual  coherence  functions  T see  Eq.  (33)  in  Topical  Report  1]. 

The  principal  parameter  is  the  size  of  the  aperture  compared  with  the 
spatial  coherence  scale. 

The  temporal  variations  of  the  signal  mainly  affect  the  sensitivity 
of  the  system  as  distinct  from  its  resolving  power.  Tn  effect,  the  target 
cross  section  fades--a  problem  that  has  been  treated  extensively  in  radar 
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theory.  ’  Similarly,  loss  of  phase  coherence  limits  the  effective 

U 

coherent  integration  time.  Finally,  loss  of  frequency  coherence  limits 
the  achievable  range  resolutions.  Such  effects  are  not  fundamentally 
different  from  their  counterparts  in  SATCOM  systems,  and  detailed  simulations 
are  ultimately  required  for  performance  evaluations. 

In  the  remainder  of  this  report,  we  consider  the  performance  degrad¬ 
ation  of  sate  1 1 i te -borne  SAR  systems.  The  SAR  method  achieves  high  res¬ 
olutions  (comparable  to  practical  optical  instruments  in  space)  by  co¬ 
herently  processing  the  radar  returns  over  periods  approaching  the  time 
a  target-element  remains  within  the  beam.  Coherent  processing  is  used 
to  compensate  for  the  phase  drift  induced  by  the  satellite  motion.  Pro¬ 
pagation  disturbances  cause  both  phase  and  amplitude  errors  that  distort 
the  synthesized  beam. 

We  have  divided  our  treatment  of  the  problem  into  three  parts.  In 
Section  II,  we  develop  a  model  that  includes  all  the  essential  elements 
of  a  SAR  system;  we  used  the  SEASAT-A  SAR  as  a  model.  The  satellite 
was  operational  only  for  a  little  more  than  three  months  commencing  in 
July  1978,  but  it  produced  a  large  data  base  for  demonstrating  propagation 
effects.  In  Section  III,  we  used  the  SATCOM  model  with  typical  iono¬ 
spheric  parameters  to  determine  the  likelihood  of  significant  SAR  per¬ 
formance  degradation  due  to  natural  disturbances. 

Our  analysis  shows  that  only  the  most  intense  high-1  at i tude  distur¬ 
bances  are  likely  to  produce  a  detectable  SAR  image  degradation.  More 
severe  disturbances  would  be  expected  near  the  geomagnetic  equator.  In 
co 1 1 abora t ion  with  staff  members  from  Research  Development  Associates 
(RDA)  who  were  working  with  SEASAT  data  for  other  reasons,  we  performed 
a  search  ot  the  SEASAT-B  library  of  optically  processed  data.  One  pass 
that  showed  a  conspicuous  image  degradation  coincided,  fortuitously, 
with  a  set  of  incohe ren t-sca t ter  radar  measurements .  The  radar  data 
showed  a  prominent  F-layer  enhancement  with  the  characteristics  that  have 
been  consistently  associated  with  enhanced  sc  in t i 1 1  a t i on . 

In  Section  V  we  discuss  the  effects  of  more  severe  propagation 
disturbances  and  future  effort  to  identify  and  mitigate  the  effects  of 
propagation  disturbances. 
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II  A  MODEL  FOR  EVALUATING  PROPAGATION  EFFECTS  IN  SPACE-BASED  SAR 
SYSTEMS 


The  principal  of  SAR  processing  is  simple.  A  linear  array  of 
discrete  elements  can  be  focused  by  adjusting  the  phase  of  the  individual 
elements  so  that  the  returns  from  a  signal  emanating  from  the  focal  point 
have  the  same  phase.  Alternatively,  the  complex  signal  from  a  single 
element  can  be  recorded  as  it  moves  along  the  axis  of  the  leer  array, 

and  the  appropriate  phase  adjustments  can  be  made  ex  post  to.  The 

beamwidth  of  a  fixed  antenna  element  imposes  a  limitation  the  maximum 
aperture  that  can  be  synthesized. 

For  example,  let  the  aperture  width  of  the  antenna  element  be  L  . 

E 

The  maximum  element  separation  that  can  contribute  to  a  focus  at  R  is 


The  spatial  resolution  of  the  synthesized  array  is 


Ay’  "  2L  ~  2 


Because  the  SAR  antenna  is  used  for  both  transmission  and  reception,  the 
effective  aperture  is  halved;  hence,  the  factor  of  2  in  Eq.  (2). 

For  the  SEASAT-A  satellite,  a  phased  array  with  an  effective  aperture 
of  A.  12  m  was  used;  thus,  the  intrinsic  resolution  is  a  6  m.  The  aperture 
is  traversed  in  2  s.  In  practice,  four  4-kin  segments  are  processed  sepa¬ 
rately  and  then  averaged  to  reduce  cross  section  variations  caused  by  the 
changing  aspect  angle.  Thus,  t  21  n  is  the  azimuth  resolution  usually 
achieved  lor  digitally  processed  SAR  images. 

High  resolution  in  the  cross-track  direction  is  achieved  by  using 
pulse  compression.  The  transmitted  signal  has  the  form 
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P  ( t )  =  Ro  i  roc  t  (TJ  exPr2«iUt+  K,.t  )  ]  j  .  (  3  ) 


For  SEASAT-A, 

T  =  33.9  s 

=  0.362  MHz  / s 

f  =11 .29  MHz 

o 

Details  of  the  SEASAT-A  signal  structure  and  data  processing  can  be 
found  in  References  9  and  6.  The  corresponding  complex  signal  will  be 
denoted  by  p ( t  )  . 

To  model  the  return  signal,  we  consider  the  coordinate  system  shown 
in  Figure  2,  where  r(t)  is  the  distance  i rom  the  satellite  at  Lime  t  to 


FIGURE  2  LOCAL  COORDINATE  SYSTEM  FOR  SAR  PROCESSING 


a  fixed  point  (x,  y)  on  the  ground-  Thus, 


r  (  t )  = 


(Rsin0  +  x)  +  (y  -  vt)  +  (Rcos0)* 


2  2 

-  D  •  o  yvt  p  +(vt) 

=  R  +  xslnQ  '  —  +  2R - 


(4) 


2  2  \ 

where  p  =  (x  +  y  )  .  The  return  signal  from  a  single  point  target  at 
p  admits  the  representation 


S ( t )  =  Re  {  p  (  x ) *g ( t ,  x  )  o(p)exp[4nir(t)/A]}  ,  (5) 


where  o(p)  is  the  complex  scattering  coefficient  assumed  to  be  frequency 
independent  over  the  frequency  band  of  the  pulse.  Ionospheric  effects 
are  accommodated  by  convolving  p(x)  with  the  time-varying  impulse  function 
g  (t,x).  Alternatively,  we  can  write 


p  ( x ) *g ( t ,  x  ) 


h 


=  lp( f  )h(t ,  f 


+  f  )exp  (  2it  if  x  )d  f 


(6) 


where  h(t,  f)  is  the  time-varying  transfer  function  and 


g(t  ,x  ) 


f  +  f )exp( 2n i f x  )df 


(  7) 


In  the  SEASAT-A  data  reduction,  bursts  of  13,680  samples  over  a 
300.46-ps  interval  are  recorded  every  interpulse  period  of  Tp  =  600  us; 
thus,  1/At  =  45  MHz.  The  spectral  extent  of  the  sampled  signal  is  shown 
in  Figure  3(a).  The  data  are  processed  by  first  performing  a  4096-point 
real  Fast  Fourier  Transform,  which  generates  2048  complex  frequency  samples 
spanning  f  =  0  to  f  =  22.5  MHz.  To  generate  the  complex  signal  ,  a  cyclic 
shift  of  the  Fourier  components  is  performed  as  illustrated  in  Figure  3(b). 

The  complex  signal  is  then  multiplied  by  p  (nAf).  Tf  the  t i mo -bandw id th 
product  is  largo,  the  inverse  Fourier  transform  of  |p|2  is  well  approx¬ 
imated  hy 


p  (t) 
c 


K  l'2s  in  (  n  K  Tt ) 
r  r 


nK  Tt 
r 


(8) 
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EXTRACT  WITH  REAL  EFT 


FIGURE  3  SIGNAL  SPECTRUM  FOR  RANGE  COMPRESSION 


In  practice,  the  radial  component  of  the  satellite  motion  causes  a 
Doppler  shift  that  manifests  itself  as  an  uncertainty  in  f  ,  and  a  time- 
consuming  iterative  procedure  must  be  used.  For  our  purposes  here,  how¬ 
ever,  this  subtlety  can  be  ignored. 

We  note,  however,  that  over  the  300-ps  data  interval,  changes  in 
r(t)  are  negligible,  as  are  changes  in  g(t,x)  with  t.  Thus,  the  pulse 
compression  operation  affects  only  p(t)  in  Eq .  (5);  moreover,  At  = 
(K^T)-^.  In  range  delay  units,  we  can  write 


(.!,*> 


■j^JAR  -  xsinQ  -  n£)*g(£Tp,  jAR  -xsin©  -  n^) 


expCAiti-^— )  cs(p)dp  , 


(9) 


where  Tp  is  the  interpulse  period;  then. 


vTp£ 

n(!  =  -y  —  +  ^ 


2  2 
p  +  ( vT  i  ) 


(10) 


and  we  have  integrated  over  all  scatterers. 

The  fact  that  appears  in  the  argument  of  p^fx)  greatly  complicates 
the  azimuth  compression  operation.  The  range  cell  for  which  the  azimuth 
compensation  n^  is  to  be  applied  changes  with  £.  To  make  this  explicit, 
we  define 


j  U )  =  l(xsin@  -n^l/ARjj  ,  (11) 

where  [  .  ]  denotes  the  nearest  integer  value.  The  processed  signal 

for  a  target  area  centered  at  (x',  y 1  )  is 


Np-1 

d(p’)  =  ^  sJj'U),  £  ]exp( -Ait  iq  '  ^  /  X  )  .  (12) 

£=0 


Because  of  the  phase 
2 

interest  is  | a( p ' ) | 


fluctuations  of  the  target  area,  the  quantity  of 
,  where  the  overbar  denotes  an  average. 
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One  usually  assumes  that  the  target-scattering  function  is  a 
random  process  that  is  uncorrelated  from  one  resolution  all  to  the  next. 
We  define  the  ambiguity  function  as 


Np-  l 


2-0 


p^*g(2Tp,  j'(2)AR  -  xsinQ  -  n£) 


x  exp{4ni(n^  -  n'^)/X} 


(13) 

Using  the  assumption  of  uncorrelated  resolution  cells,  it  follows  that 

(14) 


<|o(p')|2>  =  J'J' <  |  A(p  ,  p  '  )  |  2><  |  o(p  )  |  2>dp 


The  ionospheric  disturbance  potentially  affects  both  the  range,  x,  and 
azimuth,  y,  coordinates;  however,  for  naturally  occurring  ionospheric 
disturbances,  the  range  distortion  is  negligible,  as  we  shall  show  in  the 
next  section.  It  should  also  be  noted  that  insofar  as  evaluating  the 
effects  of  propagation  disturbances  are  concerned,  there  is  no  loss  of 
generality  in  assuming  that  the  compressed  pulse,  pc ,  was  actually 
transmitted . 
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Ill  IONOSPHERIC  EFFECTS 


The  mathematics  of  wave  propagation  in  randomly  irregular  media 

7  8 

have  been  treated  extensively  in  published  papers,  ’  as  have  its 

9 

applications  to  systems  effects.  Here  we  shall  only  review  these  re¬ 
sults  to  justify  our  choice  parameters  used  in  evaluating  the  effects  of 
naturally  occurring  ionospheric  disturbances.  Scintillation  is  essentially 
an  interference  phenomenon  as  is  illustrated  schematically  in  Figure  4. 


STRUCTURED 

REGION 


I 

P 


FIGURE  4  SCHEMATIC  REPRESENTATION  OF  SCINTILLATION  DEVELOPMENT 


The  spatial  wavenumber  spectrum  of  the  complex  field  u  (p,z) 
characterizes  the  angular  distribution  of  the  scattered  waves.  Formally 
the  wavenumber  spectrum  is  the  Fourier  transform  of  the  mutual  coherence 


function  which  has  the  particularly  simple  mathematical  form 


<u (p , z  )u*(p '  , z )>  =  exp{ -iDfAp  ) } 


(15) 


where  D(Ap)  is  the  phase  structure  function.  The  in-situ  irregularities 
are  characterized  by  a  three-dimensional  power  spectral  density  function 
of  the  form 

-(p+1 ) 


$  ( q  )  =  C^q 


(16) 


where 


C  =  8n3/2q  P'2  - - T 

s  Mo  „/p-2\ 


■(¥) 


J  P-2| 
2  / 


<an2> 

e 


(17) 


relates  the  turbulence  parameter  to  the  electron  density  variance 

via  an  outer  scale  wavenumber  q^ .  The  corresponding  phase  variance 

for  a  path  of  length  l  is 

P 


(¥) 


<6^2  >  =  2/iur  ^X^  — T~o —  q  G<AN^>]  , 
/p-2v  o  p  c 


■(¥) 


(18) 


where  r  is  the  classical  electron  radius  and  G  is  a  geometric  factor 
e 

that  accounts  for  irregularity  anisotropy.  The  square  brackets  denote 
an  average  va  1  vie  ewer  the  propagation  path. 

In  terms  of  these  model  parameters,  it  can  be  shown  that 


D(y)  =  r2X2«pCsC(p)|y|max(2,  P_1)  2<p<4  , 


(19) 


where 


C(P)  =  -¥•  f 0 .25  -2.25  p  +  5.5] 


(20) 


We  define  ?  so  that 

o 


D(e  )  =  2 

o 


(21  ) 
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The  corresponding  temporal  coherence  time  is  given  as 


*”o  v 


If  we  let  p  =  3,  Eq.  (22)  is  equivalent  to  Eq .  (2)  in  Wittwer.  In 

that  case,  Z  a  f . 

o 

The  frequency  coherence  can  be  estimated  by  noting  that  the  average 

angular  deviation  is  proportional  to  (£  k)  \  which  we  can  equate  with 

°  <502 

g 

60  in  Figure  4.  The  path  difference  for  small  60  is  z— .  The  corres- 

S  S  t- 


ponding  delay  spread  is  given  as 


,  2,2 

ck  Z 

o 


The  corresponding  frequency  spread  f^  is  defined  as 


f  =  J- 
o  At 


2  2 

ck  £ 

o 


Again,  if  we  let  p  =  3,  Eq.  (24)  is  equivalent  to  Eq .  (3)  in  Wittwer 
for  which  f  «  f  ^ .  In  any  case,  the  coherence  bandwidth  changes  much 
more  rapidly  with  frequency  than  does  the  spatial  coherence,  Z  . 

The  general  form  of  frequency  correlation  function  does  not  admit 

9 

simple  mathematical  representation.  Wittwer  has  derived  a  representation 
valid  for  p  =  3.  An  integral  representation  based  on  the  phase  screen 
model  was  used  to  interpret  data  from  the  Wideband  satellite. ^  The 
main  differences  between  these  models  lie  in  the  approximations  that 
have  been  used  and  the  parameter  ranges  that  have  been  emphasized. 


The  SATCOM  model  for  nuclear  effects  assumes  a  q  power  law  to  an 

outer  scale  wavenumber  corresponding  to  'v  10  km.  For  the  ionosphere, 

11  12 

a  two-component  power  law  seems  most  representative.  ’  For  inter¬ 
mediate-scale  structures  between  10  km  and  500  m,  a  shallowly  sloped 
power  law  with  2  <  p  <  3  applies.  For  sizes  smaller  than  500  m,  a  much 
steeper  slope,  3  <  p  <  5,  applies. 
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A  detailed  analysis  of  the  scintillation  effects  of  a  two-component 

13 

power  law  has  been  performed.  The  resultant  scintillation  structure 
depends  critically  on  the  Fresnel  radius  r^..  If  r^  is  smaller  than  the 
scale  of  the  break  frequency,  the  scintillation  structure  is  essentially 
the  same  as  a  single  component  power  law  with  a  steep  slope  (p  >  3);  r^ 
is  greater  than  the  scale  of  the  break  frequency  but  less  than  the  outer 
scale,  the  break  acts  much  like  an  outer  scale  cutoff.  Finally,  as  the 
Fresnel  radius  increases  toward  the  outer  scale,  the  low-frequency  por¬ 
tion  of  the  spectrum,  where  most  of  the  spectral  intensity  lies,  dominates 
the  signal  moments;  however,  the  high-frequency  components  of  the  scin¬ 
tillation  intensity  spectrum  follows  the  steeply  sloped  portion  of  the 
integrated  phase  spectrum. 


At  L-band  and  higher  frequencies,  typical  Fresnel  radii  are  smaller 
than  the  scale  of  the  break  frequency  q^  .  Thus,  the  scintillation 
structure  has  the  characteristics  of  a  steeply  sloped  power  law  environ¬ 
ment.  To  extend  Eq .  (16)  beyond  q^  ,  we  can  use 


®  (q  ) 


-P« 


C  q/VV 
s  b 


%  <  q  i  qb 


q  >  q. 


(25) 


It  foil ows  tha t 

C'  =  C  q, (ph'P2}  (26) 

s  s  ’b 

is  the  effective  turbulent  strength  for  the  transition  scale  structure. 

In  Figure  5  we  have  plotted  t  and  f  with  the  S.  scintillation 

o  o  4 

index  using  both  the  intermediate-scale  (p  )  and  the  transitional  scale 

e 

(p  )  power  law  components.  A  break  frequency  of  500  m  was  used,  which 

1  12 
exceeds  the  median  value  of  750  m  reported  by  Basu  et  al .  Thus,  the 

curves  in  Figure  5  should  bracket  the  actual  coherence  scales.  In 

general,  the  effect  of  the  steeply  sloped  transition  scale  is  to  increase 

both  the  coherence  bandwidth  and  the  spatial  coherence.  In  all  cases, 

however,  the  coherence  bandwidth  remains  well  beyond  the  frequency  band 
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of  the  chirp  pulse.  Under  such  conditions,  the  ionosphere  causes  only 

a  pulse-to-pulse  variation,  and  p  -"g  in  Eq .  (13)  can  be  replaced  by  the 

c 

simple  product  p^  (x  )h(j?T  ,  f  ) ,  which  follows  from  Eq .  (7). 

The  spatial  coherence,  however,  can  be  reduced  below  the  limits  of 
space-based  SAR  apertures.  To  evaluate  the  effect  of  propagation  dis¬ 
turbances  in  such  cases,  we  have  simulated  the  ambiguity  function  as 
defined  by  Eq .  (13),  using  the  parameters  listed  in  Table  1.  For  simpli¬ 
city,  we  have  ignored  the  range  migration  effect  so  that  the  propagation- 
induced  distortion  is  confined  entirely  to  the  azimuth  coordinate. 
Single-phase  screen  simulations  based  on  the  two-component  power  law 

were  used  to  generate  realizations  of  h(£T  f  ). 

K  C 


Table  1 

SEASAT-A  PARAMETERS 


Symbo 1 

Name 

Value 

f 

F  tequency 

1275  MHz 

T 

p 

[nterpulse  period 

(19(54)-’  s 

V 

Veloc i ty 

8.78  km / s 

N 

p 

Number  of  pulses 

2(5(57 

L 

Synthesized  aperture  size 

1  f>  km 

R 

Slant  range 

830  km 

0 

Beam  offset 

20° 

Ax 

Range  step 

a.  8  m  (not  used 
in  simulation) 

Ay 

Azimuth  step 

vT  =  6  m 

P 

Representative  results  for  progressively  increasing  C  values  are 
shown  in  Figures  6  through  8.  The  effect  of  the  propagation  disturbances 
is  to  lower  the  peak  of  the  main  lobe  and  raise  the  sidelobe  level. 
Broadening  of  the  main  lobe,  which  one  would  associate  with  loss  of 
resolution,  is  not  significant.  Thus  the  propagation  effects  act  to 
reduce  contrast  rather  than  to  defocus  the  SAR  image.  This  is  an  importan 
distinction  for  processing  algorithms  and  mitigation  techniques. 
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FIGURE  7  SYNTHESIZED  ANTENNA  BEAM  FOR  PROPAGATION  DISTURBANCE  WITH  C, 


Other  analyses  of  propagation  disturbances  in  SARs  have  been  made,*'4 
these  analyses  confined  the  errors  entirely  to  the  phase  of  the  signal. 

A  least-squares  estimate  of  the  linear  and  quadratic  phase  terms  was 
obtained,  together  with  rms  value  of  the  residual  called  the  integrated 
side  lobe  ratio.  Consistent  with  our  own  results,  the  integrated  side- 
lobe  ratio  was  found  to  lie  the  most  serious  concern  tor  SAR  imago 
degrada  t ion. 
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rv  SEAS AT- A  DATA 


A  large  number  of  high-latitude  SAR  images  were  recorded  by  the 
SEASAT-A  satellite  when  it  was  operational  in  1978.  A  systematic  searcli 
of  the  SEASAT  library  of  passes  was  performed  to  identify  pass  segments 
where  propagation  disturbances  were  likely.  In  the  initial  sorting,  we 
accepted  all  pass  segments  that  fell  within  the  auroral  zone  or  polar 
cap  within  four  hours  of  local  midnight.  The  auroral-zone  passes  were 
then  further  divided  into  nighttime  Alaska  sector  and  nighttime  non- 
Alaska  sector  passes. 

The  SEASAT  passes  that  satisfy  these  criteria  are  summarized  in 
TabLe  2.  These  passes  collected  data  during  the  most  likely  periods  for 
auroral  activity.  We  made  no  attempt  to  further  sort  the  data  on  the 
basis  of  magnetic  activity  or  other  indicators  of  auroral  activity; 
rather,  we  decided  to  make  a  systematic  visual  survey  of  the  optically 
processed  SEASAT  images  to  identify  any  degradation  that  could  be 
attributed  to  auroral  propagation  disturbances. 

Because  several  factors  that  are  difficult  to  control  affect  the 
quality  of  the  optical  images,  we  sought  distinct  changes  within  a  single 
image  that  could  not  he  immediately  associated  with  actual  changes  in  the 
terrain  and  other  factors.  We  paid  careful  attention  to  operator  logs 
and  comments  by  personnel  familiar  with  the  detailed  processing  techniques. 

Dr.  Victor  Gonzalez  of  SRT  and  Dr.  Dean  Liskow  of  RAD  paid  two  visits 
to  the  let  Propulsion  Laboratory  (IPL).  About  (SO  of  the  optically  pro¬ 
cessed  SEASAT  passes  from  the  candidate  subset  were  examined.  The 
optically  reduced  data  are  presented  in  the  form  of  four  strips  of 
(negative  and  positive  prints)  film,  each  several  feet  long.  The  quality 
of  the  data  is  good  most  of  the  time,  and  the  images  are  clear. 

The  resolution  in  the  data  studied  was  of  the  order  of  10  to  20  km 
per  inch.  Thus,  the  resolution  did  not  allow  examination  of  small  de¬ 
tails  of  less  than  about  100  m.  Rivers,  mountain  ranges,  and  coast 
lines  provided  good  edges  for  comparison  during  different  passes  because 
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Table  2 


CONDI  DATE  PASSES  FOR  AURORAL  PROPAGATION  DISTURBANCES 


(a  ) 

Alaska  Auroral  Zone 
(primary  candidate) 


(b> 

Polar  Cap 

(supplementary  passes) 


( c  ) 

Non-Alaska  Auroral  Zone 
(Last  Resort) 


ID  « 

Start 

IUT) 

Stop 

(UT) 

001 

mm 

n 

002 

IpppB 

■rH 

004 

H 

1338 

013 

1200 

1206 

028 

1246 

1232 

03  7 

1233 

1239 

040 

1401 

140  7 

04b 

1331 

1  337 

032 

141  1 

1413 

099 

1  329 

1333 

106 

1437 

1443 

1  46 

1437 

1  303 

13] 

1428 

1434 

163 

143  3 

1441 

1  74 

1439 

1  44  3 

188 

1  634 

1  700 

202 

1  320 

1  624, 

303 

1026 

1031? 

3  30 

1039 

1  04  3  “ 

333 

1032 

1  03  3 

419 

1  336 

1  346 

-+2  7 

1  129 

1  1  34 

■4-42 

1  349 

1  339 

44 '3 

1  703 

1  713 

4  30 

1142 

1  147 

46  5 

1402 

141  1 

468 

1  71  7 

1  72b 

4  74 

1  133 

1  202 

Mer  id i an 


ID  « 

Stop 

(UT) 

086 

092 

09  7 
104 

1  19 
124 
128 
132 
137 

142 

143 
149 

161 

169 

1  73 

181 

193 

199 

2  76 
280 

00 

00 

283 

CO 

< 

on 

< 

283 

cc 

288 

296 

r) 

u* 

ro 

u- 

304 

ID 

U2 

315 

r~' 

ZJ 

a 

322 

cj 

o 

340 

z 

1 — ' 

349 

357 

366 

3  74 
381 
389 
397 
404 
412 
428 
436 
451 
459 

4  76 

ID  * 

S  tart 
(UT) 

Stop 

(UT) 

016 

1130 

1136 

025 

1000 

1006 

032 

103  7 

1043 

038 

1432 

1438 

044 

1014 

1021 

045 

1152 

1157 

048 

0944 

0951 

089 

2125 

2131 

107 

2132 

2137 

121 

2139 

2145 

129 

1343 

1349 

143 

1351 

1356 

153 

2231 

223  7 

165 

2237 

2243 

1  76 

2134 

2140 

218 

2340 

2342 

233 

2353 

2354 

248 

0005 

000  7 

264 

001  7 

0020 

2  74 

0030 

0033 

308 

2353 

2357 

319 

0108 

0111 

344 

0008 

0011 

344 

0121 

0124 

3  52 

0230 

0232 

361 

0021 

0023 

3  75 

1626 

1632 

377 

0103 

0107 

398 

1639 

1  64  5 

400 

0116 

0120 

422 

1652 

1  701 

439 

0031 

0037 

462 

0044 

004  5 

480 

0124 

0130 
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Che  data  were  frequently  taken  over  the  same  area.  In  some  instances, 
however,  the  configuration  of  the  ice  in  rivers  changed  enough  to  make 
the  recognition  of  details  difficult. 

Examination  of  the  pictures  showed  that  they  contained  times  (or 
areas)  in  which  the  data  were  deficient  to  various  degrees,  and  sometimes 
were  missing  completely.  Discussion  with  JPL  staff  familiar  with  the 
acquisition  of  data  disclosed  that  several  reasons  could  account  for  data 
degradation.  For  example,  dynamic  behavior  of  the  satellite  could  be 
one  reason.  When  the  data  are  transmitted  to  the  ground-monitoring 
stations,  the  ground  receivers  have  antennas  with  different  gains  and 
receivers  with  different  characteristics  that  add  "noise"  to  the  analog 
data  received.  The  details  of  the  optical  processor  also  have  adjustments 
(such  as  focusing)  that  also  may  occasionally  add  noise  to  the  data. 

Nonetheless,  we  isolated  a  pass  that  did  contain  the  expected 
characteristics  of  a  propagation  disturbance: 

ID  No.  330 

Revolution  No.  1236 

Date  21  Sept  1979 

UT  1061:50  to  1062:10. 

We  also  located  an  unperturbed  pass  that  covered  the  same  area: 

ID  No.  355 

Revolution  No.  1279 

Date  26  Sept  1979 

UT  1056:33  to  1056:33. 


The  characteristics  of  the  degradation  arc  as  follows: 
•  The  degraded  area  is  localized. 


The  degraded  area  does  not  include  all  ranges. 

The  boundaries  of  the  degraded  area  are  not  perpendicular 
to  the  direction  of  satellite  motion. 
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Figure  9  shows  the  image  swatii  of  the  two  SEASAT-A  passes  superimposed 
on  a  map  of  Alaska  and  northern  Canada;  the  shaded  area  is  the  region  of 
interest.  Figure  10  shows  a  topographical  map  of  the  image  area,  which 
lies  northeast  of  Fort  Yukon,  Alaska. 

Figure  11  shows  the  digitally  processed  normal  SAR  image.  The 
Porcupine  and  the  Coleen  Rivers,  as  well  as  the  nearby  lakes,  correlate 
very  well  with  the  topographical  map.  The  striking  map-like  character 
of  SAR  radar  images  is  well  known.  Figure  12  Is  a  radar  image  of  the 
same  area  but  recorded  several  days  earlier.  This  image  shows  a  broad 
wedge-shaped  region  where  there  is  a  distinct  loss  of  contrast  but  no 
perceptible  loss  of  resolution.  This  type  of  degradation  is  consistent 
with  the  analysis  presented  above  in  Section.  III. 

Through  an  extremely  fortuitous  set  of  circumstances,  the  Chatanika 
radar,  which  is  located  near  Fairbanks,  Alaska,  was  operating  during 
the  period  the  SAR  image  shown  in  Figure  12  was  recorded.  The  radar  was 
performing  a  routine  24-hour  run  in  a  mode  that  included  elevation  scans 
as  well  as  multiposition  measurements  for  drifts.  From  the  elevation 
scan  data,  it  is  possible  to  map  the  distribution  of  ionization  in  the 
magnetic  meridian.  Such  elevation  scans  have  been  used  extensively  to 
study  unstable  F-region  "blobs";  indeed,  they  provided  the  first  demon- 

1  5 

stration  of  the  association  of  F-rcgion  structure  and  scintillation. 

Figure  13  shows  the  track  of  the  subsatellite  point  and  a  tracu  dis- 
P'  aced  20°  from  the  orbit  plane  as  viewed  from  the  satellite.  The  points 
on  the  subsatellite  track  are  20  s  apart.  The  numbers  on  the  20°  track 
are  the  geometrical  scintillation  enhancement  factors  for  10:10:1  sheets.' 
Tt  is  interesting  to  note  that  this  pass  is  not  favorably  oriented  for  a 
strong  geometrical  enhancement  in  the  image  region. 

The  dark  bar  locates  the  meridian  plane  of  the  radar  scan.  The 
area  of  the  SAR  image  is  outlined,  and  the  region  of  the  degradation  is 
shaded.  The  plane  of  the  meridian  scan  intercepts  the  image  region; 
however,  the  F-region  penetration  point  of  the  structure  region  lies  to 
the  east  of  the  ground  image  area.  It  is  significant  that  the  degradation 
was  observed  as  an  east-west  protrusion  into  the  image  region. 
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FIGURE  12  DIGITALLY  PROCESSED  IMAGE  FOR  REVOLUTION  1236 
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The  dark  bar  locates  the  meridian  plane  of  the  radar  scan.  The 
area  of  the  SaR  image  is  outlined,  and  the  region  of  the  degradation  is 
shaded.  The  plane  of  the  meridian  scan  intercepts  the  image  region; 
however,  the  F-region  penetration  point  of  the  structure  region  lies  to 
the  east  of  the  ground  image  area.  It  is  significant  that  the  degradation 
was  observed  as  an  east-west  protrusion  into  the  image  region. 
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The  satellite  track  identifier  closest  to  the  image  region  occurred 
at  10:30.  Two  of  the  12-minute  elevation  scans  bracketed  this  period. 

The  reconstructed  meridional  density  contour  maps  are  shown  in  Figures 
14  and  15. 

A  highly  structured  enhanced  F-layer  is  a  prominent  feature  in  both 
maps.  From  recent  studies  of  these  F-region  structures,  we  know  that 
they  are  highly  unstable;  moreover,  although  the  large-scale  structure 
maps  along  the  L-shells,  there  is  a  considerable  amount  of  small-scale 
east-west  structure  as  predicted  by  Keskinen  and  Ossakow.^  Thus,  we 
cannot  correlate  the  individual  blob  structures  with  the  image  degradation. 
However,  we  know  that  scintillation-producing  irregularities,  of  the  type 
we  have  been  studying  intensively  for  several  years  now,  were  almost 
certainly  present  along  the  propagation  path  to  the  degraded  image. 

We  have  emphasized  the  data  from  SEASAT-A  revolutions  1236  and  1279 
because  of  the  incoherent-scatter  radar  data  that  allowed  us  to  identify 
the  likely  source  region  of  the  disturbance.  Also,  in  our  initial  search 
of  the  candidate  passes  listed  in  Table  2,  we  were  looking  for  image 
degradation  of  an  unknown  character.  The  tendency  in  the  initial  search 
was  to  look  for  loss  of  resolution  rather  than  loss  of  contrast.  The 
optically  processed  images  that  were  used  for  the  initial  search,  more¬ 
over,  were  processed  in  four  strips,  which  often  showed  intensity  gradients 
across  the  image. 

Thus,  it  is  reasonable  to  assume  that  3  propagation-induced  contrast 
reduction  in  h i gh-la t i tude  SAR  images  is  present  in  many  more  of  the 
candidate  passes  but  was  overlooked  in  the  initial  search.  Indeed,  Dr. 
Bryan  Honeycutt  of  1PL  had  noted  a  contrast  reduction  in  a  SEASAT-A  image 
of  Banks  Tsland  during  a  period  of  high-geomagnetic  activity.  He  attri¬ 
buted  the  effect  to  an  aurora  1 -ac t  iv i ty- induced  propagation  disturbance. 
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INVARIANT  LATITUDE 


V  CONCLUSIONS  AND  RECOMMENDATIONS 


In  this  report  we  have  analyzed  the  likely  effects  of  propagation 
disturbances  on  satellite-borne  SAR  imagery.  The  work  is  an  application 
ot  the  general  space-hased  radar  model  developed  during  the  first  phase 
ol  this  contract.  The  principal  effect  of  propagation  disturbances  is 
a  degradation  of  the  azimuthal  component  of  the  SAR  ambiguity  function. 

In  a  very  severely  disturbed  environment,  coherence  bandwidth  loss  could 
also  degrade  the  range  resolution,  although  at  L-band  and  higher  frequen¬ 
cies  naturally  occurring  disturbances  will  not  cause  such  effects. 

Simulations  of  the  azimuthal  degradation  show  that  an  elevated 
sidelobe  level  is  the  dominant  effect.  Tn  an  actual  SAR  image,  the  ele¬ 
vated  sidelobes  evidently  reduce  contrast.  We  showed  a  very  good  example 
of  this  effect  in  an  image  obtained  from  a  SEASAT-A  pass  over  Alaska  with 
correlative  incoherent -scatter  radar  data  to  identify  the  disturbance. 

More  severe  disturbances,  of  the  type  commonly  observed  near  the 
geomagnetic  equator,  would  likely  cause  a  much  more  severe  contrast  re¬ 
duction,  although  the  detailed  characteristics  of  such  disturbances  should 
be  investigated  further. 

A  potentially  important  factor  in  the  SKASAT-A  image  we  analyzed 
is  that  the  propagation  disturbance  occurred  in  a  geometry  for  which 
the  enhancement  of  the  disturbance  and  to  the  known  sheet-like  anisotropy 
of  the  irregularities  was  not  significant.  Thus,  f or  a  more  nearly 
meriodinal  alignment  of  the  satellite,  the  same  structure  could  produce 
a  more  severe  degradation.  More  elaborate  simulations  than  those  we 
performed  in  Section  III  should  be  initiated  to  identify  the  effects  of 
the  geometric  dependence  of  propagation  disturbances. 

Such  simulations  a^so  would  provide  a  data  source  for  testing  more 
elaborate  mitigation  schemes.  As  shown  in  Section  TI,  the  simulations 
need  not  reproduce  in  complete  detail  the  SAR  data  recording  and  pro¬ 
cessing  procedures.  For  example,  the  pulse  compression  operation  can 
he  eliminated  by  assuming  that  the  compressed  pulse  was  actually  trans¬ 
mitted.  There  is  no  loss  of  generality  insofar  as  propagation  effects  are 
concerned  . 


Our  preliminary  analysis  suggests  that  the  degradation  depends  on 
the  spatial  coherence  scale  l  which  varies  essentially  linearly  with 
frequency  and  that  naturally  occurring  irregularities  in  the  auroral  zone 
and,  par t  icu 1 ar  1  y ,  near  the  geomagnetic  equator  can  degrade  SAR  images. 

On  the  other  hand,  a  thorough  search  of  the  SEASAT-A  library  of  high 
latitude  SAR  data  did  not  reveal  many  examples  of  the  expected  effects. 

One  explanation  is  that  the  disturbed  regions  are  patchy  such  that 
the  average  occurrence  of  scintillation  based  on  long-Ler"  morphological 
data  exaggerates  the  likelihood  of  intense  irregularities  occurring  at  a 
single  point.  Alternatively,  the  focusing  process  used  in  the  SAR  image 
reconstruction  may  compensate  for  a  large  amount  (if  the  ionosphere- i nduced 
phase  distortion. 

To  investigate  this  possibility,  we  propose  more  detailed  simulations 
that  can  accommodate  two-dimensional  images.  This  will  also  allow  us  to 
evaluate  the  various  measures  of  image  degradation  that  are  being  used 
and  determine  the  improvements  that  can  be  achieved  by  increasing  fre¬ 
quencies  and/or  using  image  enhancement  techniques. 
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Pac i f ic-Sierra  Rsch  Corp 
ATTN:  E .  Field,  Jr 
ATTN:  F.  Thomas 
ATTN:  H.  Brode,  Chairman  SAGE 

Pennsylvania  State  University 

ATTN:  Ionospheric  Rsch  Lab 

Photometries,  Inc 

ATTN:  1.  Kofsky 

Physical  dynamics,  Inc 


ATTN 

E.  Fremouw 

Physical  Rsch,  Inc 

ATTN 

R.  Deliberis 

ATTN 

T.  Stephens 

R&D  Associates 

ATTN 

R.  Turco 

ATTN 

W.  Karzas 

ATTN 

H.  Ory 

ATTN 

C.  Greifinger 

ATTN 

M.  Gantsweg 

ATTN 

F.  Gilmore 

ATTN 

W.  Wright 

ATTN 

P.  Haas 

RSD  Assoc iates 

ATTN 

B.  Yoon 

Rand  Corp 

ATTN 

E.  Bedrozian 

ATTN 

C.  Crain 

ATTN 

P.  Oavis 

Rand  Corp 

ATTN 

B.  Bennett 

Riverside 

Rsch  Institute 

ATTN 

V.  Trapani 

Rockwel 1 

nternational  Corp 

ATTN 

R.  Buckner 

Rockwel 1 

nternational  Corp 

ATTN 

S .  Qu  i  1  i  c i 

DEPARTMENT  OF  DEFENSE  CONTRAaORS  (Continued; 

Science  Applications,  Inc 
ATTN:  M.  Cross 

Science  Applications,  Inc 
ATTN:  J.  Cockayne 


SRI  International 


ATTN 

R.  Tsunoda 

ATTN 

J.  Vickrey 

ATTN 

W.  Chesnut 

ATTN 

R.  Leadabrand 

ATTN 

R.  Livingston 

ATTN 

D.  McDaniels 

ATTN 

M.  Baron 

ATTN 

G.  Price 

ATTN 

D.  Neil  son 

ATTN 

A.  Burns 

ATTN 

W.  Jaye 

ATTN 

J.  Petrickes 

ATTN 

G .  Smith 

4  cy  ATTN 

C.  Rino 

4  cy  ATTN 

V.  Gonzales 

Stewart  Radiance  Lab 

ATTN:  J.  Ulwick 


Technology  International  Corp 
ATTN:  W.  Boquist 

Toyon  Rsch  Corp 

ATTN:  J.  Ise 
ATTN:  J.  Garbarino 

Tri-Corn,  Inc 

ATTN;  D.  Murray 

TRW  Electronics  A  Defense  Sector 
ATTN:  R.  Plebuch 
ATTN:  G.  Kirchner 

Utah  State  University 

Attention  Sec  Control  Ofc  for 

ATTN:  K.  Baker,  Dir  Atmos  &  Space  Sci 
ATTN:  L.  Jensen,  Elec  Eng  Dept 
ATTN:  D.  Burt 
ATTN:  A.  Steed 

VisiDyne,  Inc 

ATTN:  C.  Humphrey 
ATTN:  C.  Shepard 
ATTN:  W.  Reidy 
ATTN:  J.  Carpenter 

Science  Applications,  Inc 


ATTN: 

L. 

Linson 

ATTN: 

D. 

Haml  in 

ATTN: 

C. 

Smith 

ATTN: 

E. 

Straker 
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